A real-space description of the glass transition based on heterogeneities and entropy 
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An alternative scenario for the glass transition based on the cooperative nature of nucleation 
processes and the role of entropic effects is presented. The new ingredient is to relate the dissipation 
during the relaxation process to the release of strain energy driven by the nucleation of progressively 
larger cooperative spatial regions. Using an equiprobability hypothesis for the transition between 
different metastable configurations, we obtain a relation between the free energy dissipation rate 
and the size of the largest cooperative regions. This leads to a new phenomenological relation 
between the largest relaxation time in the supercooled liquid phase and the effective temperature. 
This differs from the classical Adam-Gibbs relation in that predicts no divergence of the primary 
relaxation time at the Kauzmann temperature but a crossover from fragile to strong behavior. 
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Introduction. After many decades of efforts, a complete 
and unifying description of the glass transition problem 
is still not available. The standard approaches to the 
glass transition, have been largely based on hydrodynamic 
(such as MCTlI) or thermodynamic entropic theories like 
that proposed by Adam, Gibbs and Di Marzio (hereafter 
referred to as AGM) nearly 50 years ago@. Surprisingly, 
still nowadays the ideal AGM theory remains not ac- 
cepted nor disproved. A salient outcome of the AGM 
theory is the prediction of the existence of a second-order 
phase transition driven by the collapse of the configura- 
tional entropy (also called complexity). It is known that 
an unambiguous definition of the complexity is possible 
only in the framework of mean-field theories where phase 
space splits into ergodic components of infinite lifetime. 
Recent approaches to the glass transition problem from 
the perspective of disordered systems have, to a large ex- 
tent, validated this mean-field scenario!!. Nevertheless, a 
complete understanding of the glass transition must go 
beyond mean field by including nucleation processes into 
any valuable theory. This poses the question about how 
the mean-field picture for the glass transition is modi- 
fied in the presence of real-space effects. In this paper 
we propose an alternative scenario for the glass transi- 
tion where spatial effects such as heterogeneities play a 
crucial role in the theoretical description of the glassy 
state being also the necessary ingredient to validate and 
reinforce some aspects of entropic mean-field approaches. 
We present a phenomenological real-space description of 
nucleation processes of cooperative nature in the spirit 
of the AGM theory. This should complement other ap- 
proaches such as mode-coupling theoryu. The main out- 
come of our analysis is that the standard relation be- 
tween the relaxation time and the complexity originally 
proposed in the AGM theory is modified and the pri- 
mary relaxation time does not depend directly on the 
complexity but on the so called effective temperature (to 
be defined later). This predicts a crossover from fragile to 



strong behavior and a saturation (and not a divergence) 
of the effective barrier or primary relaxation time at the 
Kauzmann temperature. 

A phenomenological description of heterogeneities in 
glasses. One of the most intriguing features in glasses is 
the existence of heterogeneous structuresu. Experimen- 
tally these manifest as some set of atoms which have a 
dynamics manifestly slow when compared to the rest. Al- 
though experiments or numerical simulations in this field 
are very recent, heterogeneities are a direct manifestation 
of the cooperative nature of nucleation processesU. The 
basic idea in our approach is to assume that nucleation 
processes take place everywhere inside the glass when 
some structures of size s are built up by a cooperative 
mechanism. To built these structures the s atoms of the 
cooperative region must coherently move to occupy cer- 
tain positions which enable that region to release some 
strain energy. Each of these moves constitutes an ele- 
mentary activated process, for instance, the exchange be- 
tween two neighboring particles. Therefore, the hetero- 
geneities observed in the experiments are transient frozen 
structures which eventually nucleate in time scales larger 
than the observational time. When a droplet of liquid nu- 
cleates it changes to a new locally disordered structure. 
The local structure of the glass is always that of a liquid 
and there is no coarsening of a given pattern whatsoeverfj. 
In our scenario the time to activate a region of size s is, 
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where r* = tq exp(B/T) is the activated time to activate 
one atom, tq being a microscopic time and B the corre- 
sponding energy barrier. Let us introduce the quantity 
n s (t) as the number of cooperative regions (we will refer 
to them as domains) of size s at time t. Experimentally 
it is well known that time correlations in the glass state 
are stretched but decay faster than any power law. Ac- 
cording to (0) this means that the distribution n s must 



1 



abruptly fall down beyond a cuttof size s* in such a way 
that long time nucleations occur with a negligible prob- 
ability. Consequently, n s (t) for s > s* must be nearly 
zero. This approach is similar in spirit to the mosaic 
theory developed by Xia and WolyneaJ who also con- 
sidered the existence of a cutoff size s* . The existence 
of this cutoff is tightly related to the cooperative char- 
acter of the dynamics itself eq. (0) and can be illustrated 
within a simple domain aggregation model. We note that 
the model we present here is oversimplified, our aim be- 
ing only to stress the relation between cooperativeness 
and the presence of a cutoff size. Let us imagine a liquid 
that is quenched to a temperature Tf where equilibration 
cannot be achieved within the time scale of the experi- 
ment. Let us think about the molecules inside the glass as 
grouped into non-overlapping domains of different sizes 
s. After nucleating, domains of a given size s destabilize 
breaking into smaller domains. In the simplest scenario 
the aggregation and breaking of domains occurs inside a 
bath of particles. We assume that, after nucleation of a 
domain of size s, domains can gain or loose one particle 
with respective probabilities g s ,h with g s + l s = p s . p s 
is the probability that a nucleation occurs which we take 
proportional to l/r(s) where r(s) is given by eq.(|f|) with 
B = tq = 1 and T — Tf. For sake of simplicity we take 
g s = gp s ,l s = IPs with g + I = 1. Consequently, the 
balance equations involve the following "chemical reac- 
tions" : T> s — » T>s—i + P with rate l s and T> s + p — > T^a+i 
with rate g s , where T> s denotes a domain of size s and p 
a single particle of the bath. The balance equations read 
(s > 2), 

dn s (t) 

— J— = l s+1 n s+1 (t) +p s _in s _i(t) -p s n s (t) . (2) 



This set of equations must be supplemented with the 
dynamical equation for the "bath" of particles which is 
deduced from mass conservation ^2 s= i sn s (t) = const. 
In this simple model there are three parameters g, I and 
(3 = 1/7/ entering the equations. By appropriately 
rescaling the time only two of them (f3,g/l) are free pa- 
rameters. Although all possible values of g, I are possible 
the interesting regime is obtained for g/l small. Phys- 
ically this means that, during the nucleation, domains 
have more probability to loose particles than to capture 
them. This is a very reasonable assumption: just before 
the nucleation takes place the domain is in an unstable 
configuration and loosing a particle seems a more prob- 
able event. We numerically solved (|J) by numerically 
integrating them using a second order Euler algorithm. 
In the left panel of figure I we show the time evolution 
for n s {t). At any time it displays a well defined time de- 
pendent cutoff value s*(t) above which the n s (t) drops 
to zero very fast. From the knowledge of n s (t), assuming 
independent exponential relaxations for the different do- 
mains, we can also obtain the two-time correlation func- 
tion, 



C(t, t + t')=^2 sn s (t') exp(-t/r(s)) (3) 

S>1 

where t denotes the time after the quenching and 
t(s) is given by expression (Q). The correlations can 
be excellently fitted by a stretched exponential law 
with a ^-dependent stretching exponent /3 S , Ct(t') — 
exp(— (t'/r^))^*). In the right panel of figure I we 
show the corresponding correlation functions as well as 
the best fits. Note that the average relaxation time r av = 
J °° dt'Ct(t') does not necessarily scale like t. A careful 
examination of the time dependence of the distribution 
n s (t), reveals that it scales like n s (t) = (1 / s*)n(s / s*) 
with t — exp(/3s*), s*(t) being the time dependent cutoff 
size. Consequently, relaxation to equilibrium is driven 
by the growth of the largest domain of size s*(t). When 
a glass is relaxing to equilibrium, the average release of 
energy to the thermal bath is done by nucleation of the 
largest domains of size s*, meaning that the advance mo- 
tion of the front of n s (t) located at s = s* is the leading 
source of energy dissipation. Nucleation processes involv- 
ing regions of size smaller than s* always occur but do 
not yield a net energy current flow to the thermal bath. 
Metastable equilibrium is reached when the cutoff size 
s* saturates to a finite value and the net energy flow 
between the glass and the bath vanishes. 




FIG. 1. n s (t) and C(t,t + t') for different times 
t — I0 14 — 10 33 (in adimensional units) for the numerical 
solution of model (^) with I = 8, g = 1 and f3 = 2.2. The pa- 
rameter r(t) and the stretching exponent are very well fitted 
by r(t) = 2.2t - 35 , p s (t) = 0.34 + 0.45t" - 06 . 



Let us now focus our attention in the energy dissi- 
pation rate in this relaxation process. This is given by 
typical free energy variation of the relaxing domains of 
size s*(t) multiplied by their number N s * (t) = Vn s * (t) = 
(V/s*)n(l), divided by their typical nucleation time t = 
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Now, the crucial difference between nucleation pro- 
cesses in glasses and other systems is the fact that nucla- 
ation in a glass occurs between metastable liquid dropso. 
The typical release of energy AF* released after nucleat- 
ing a metastable liquid drop of size s* into another drop 
does not scale with its surface (like in coarsening systems) 
nor its volume (like in standard liquid-solid nucleations) 
but is finite and independent of s* yielding, 



1 dF 
V~dt 



s*t 



(5) 



where Ap is the average free energy changea. The growth 
of s* as a function of time is s*(t) ~ T\og(t) and stops 
when n s (t = £ oq ) reaches the stationary distribution. In 
principle the value of the equilibration time t cq depends 
on the particular model under consideration. From a mi- 
croscopic theory this time can be deduced only from the 
hydrodynamic equations. We see below how we can cir- 
cumvent the hydrodynamic description by relating the 
relaxation time of the supercooled liquid to the complex- 
ity. 

The fluctuation formula. One of the main problems in 
the theory of glasses is to define the complexity S c . In 
the original AGM theory, S c was defined as the configu- 
rational part of the entropy obtained by disregarding the 
kinetic part of the energy. Here we consider a different 
dynamical definitioB||-adapted from concepts developed 
in spin-glass theoryBtl The key concept in our defini- 
tion is the concept of conformation or basin: a basin 
includes those configurations which do not release an ex- 
tensive (proportional to s) strain energy to the thermal 
bath after nucleating a region of size s. Let us consider 
a relaxing liquid quenched to Tf. We define— S c (t, F) as 
the logarithm of the total number of basinsE£l with free 
energy F which can release strain energy to the thermal 
bath after nucleating regions of size s* at time t. In other 
words, because s* is the maximum size of nucleating re- 
gions at time t, S c (t,F) counts the number of still not 
visited conformations at time t. In this scenario regions 
of size s which have already nucleated once (i.e. s < s*) 
have already released their strain energy while regions 
which have not yet nucleated (s > s*) still contain some 
strain energy. Nucleations which do not lead to a new 
release of strain energy do not yield new conformations. 

Each of these conformations at time t may 
contain different possible configurations C which 
do not contribute to the complexity but con- 
tribute to the free energy of the conformation B, 
F B (t,T f ) = -T f \og(J2 CeB exp(-E(C)/T f ) where time 
dependent expectation values of a given observ- 
able Tf) are defined as follows: Ag(t,Tf) = 



£ Ce8 A(C)exp(- 



(E(C)-F B (t.T s )) 



The liquid character of the glass phase implies that 
basins with identical free energy Fg(t,Tf) have repro- 
ducible physical properties (i.e., independent of B). Note 



that the free energy and the expectation value of any ob- 
servable A evaluated at a given conformation B depend 
on both Tf and the quenching time t. The dependence 
on t directly appears through the set of configurations 
C contained in B. In particular, for large enough times, 
basins with very high free energy only contain configu- 
rations which have been explored several times, hence 
they do not contribute to the complexity. Contrarily, 
conformations with very low free energy contain config- 
urations which have still not been reached, hence they 
contribute to the complexity. This leads to a dynamical 
complexity S c (t, F) which has a t-dependent cutoff value 
F* such that the number of conformations or basins with 
F$(t,Tf) > F* vanishes. Moreover, the complexity is a 
monotonous increasing function of the free energy since 
the number of possible conformations which decrease the 
energy N s *(t) = (V/s*)n(l) is larger when nucleating re- 
gions are smaller. In the asymptotic large t limit it takes 
the form, 

S c (t,F < F*) = S c (F,T f ) ; S c (t,F> F*) - -oo (6) 

Now we can introduce the fluctuation formula and see 
how the present scenario substantially differs from the 
classical one by AGM. At time t after the quenching, 
nucleation processes inside the glass occur between co- 
operative regions with a disordered local structure and 
no characteristic pattern grows with time. This amounts 
to say that nucleations must be entropically driven mean- 
ing that some sort of equiprobability hypothesis for vis- 
iting conformations holds. While this is an assump- 
tion, its validity must be founded on principles simi- 
lar to those which justify the equiprobability hypothe- 
sis in Boltzmann-Gibbs equilibrium theory. Therefore, 
the probability to jump from the free energy level F* to 
another conformation with free energy F is always pro- 
portional to the number of configurations with final free 
energy F, 



Wf'^f oc 



n(F*) 



exp(Sc(F,T f )-S c (F*,T f )) . (7) 



Substituting @ in (Q) and denoting SF = F - F* we 
obtain for the transition probability, 



W SF oc e^(p sS {F*)5F)6{-5F) 



(8) 



where we have defined the effective temperature 

T cS (F*) = l//3 cS (F*) with f3 cS (F*) = ( d % F / f) ) F=F , 
being the density of complexity with free energy F*. 
This formula establishes the probability of free energy 
jumps in the aging state. Fluctuations to conformations 
or basins which increase the free energy are forbidden, 
simply because they have already nucleated in the past. 
While fluctuations to very low free energy conformations 
are entropically suppressed due to the monotonous in- 
creasing property of S c (F,Tf). In the off-equilibrium 
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state the average rate variation of the free energy at time 
t is then given by, 



1 dF X!^, xW x dx _ 1 
V~dt ~ ~tj™ oo W x dx ~ ~/3 eS (F*)t 



(9) 



We can now establish the connection between equa- 
tion ([5]) and the present one. Consistency requires that 
s* = Ap/3 c s(F*), i.e. the size of the largest nucleating 
regions s* is directly proportional to the inverse effective 
temperature (3 c s(F*) evaluated at the time-dependent 
free energy F*. After equilibrating at time t eq , F* con- 
verges to -Fcq and s* has saturated to a value s*(t eq ) 
which determines the new phenomenological relation, 



T exp^- 



T 



f 
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We can nowj-compare our prediction with the standard 
AGM scenarica. In that scenario the number of different 
conformations at time t, fl t corresponds to all possible 
combinations obtained from the two possibilities (nucle- 
ated and not yet nucleated) for the largest nucleating re- 
gions n s *(t). This number is given by f2t = 2~ yielding 
for the complexity, S c (t) = log(fi t ) = V(log(2)/s*) which 
relates the size of the cooperative region to the complex- 
ity. Using eq. (|l|) and taking t = t eq this yields the famous 
Adams-Gibbs relation r a = exp(Blog(2)V/T f S c ). The 
strong assumption contained in the AGM relation is to 
suppose that, in the off-equilibrium regime during the re- 
laxation, the glass explores nucleated and not nucleated 
conformations with the same probability. But this can- 
not be true if free energy jumps are biased towards lower 
free energy conformations. Actually, according to AGM, 
the size of the cooperative region increases with time and 
diverges at the Kauzmann temperature Tk while in our 
theory it saturates to a finite value. 

One more consequence of the fluctuation formula (||) 
concerns the fluctuation-dissipation ratio (FDR) and its 
one-step characterE2rtl After quenching to Tf a pos- 
sible way to quantify violation of FDR is to measure 
the average value of any observable A after a pertur- 
bation field Ha conjugated to the observable A is ap- 
plied to the system. Due to the disordered structure 
of the cooperative regions, if the perturbation equally 
weights all conformations, the entropically driven as- 
sumption implies that states with free energy F are sam- 
pled with a probability proportional to their number. 
In the presence of a field Ka the complexity must be 
a function of three variables, S c (F,Tf,A). We can sim- 
ply obtain the average change in the expectation value 
< A{t) > after switching on the perturbation field Ha at 
time t w . In the linear response regimeEJ the Onsager re- 
gression principle implies for the transition probabilities 
Wifm-a-M = Wsf exp(S c (F, T fl A)-S c (F, T f ,A )) = 

Wsf exp( a5 g ( ^' A) - )^ . Using the relation { 9Sc qa A) )a = 
—{3 c f[(F)hA and expanding for Ha small we finally get the 
famous violation FDT expression, g^^^ = Pcs(F*)(< 



A(t)(A(t) — A(t w )) >h A =o- Note that the description of 
the violation of FDT in terms of a single time scale t is 
consequence of the asymmetric shape (i.e. Wsf>q = 0) 
of®. 

Maiv— implications of the present scenario. It can be 
provenB that f3 cS {F) < 1/T K where S C (T K ) = at T K . 
This implies an asymptotic crossover for all fragile liq- 
uids to strong behavior. A strong glass is a fragile one 
which has exhausted all its complexity and the effective 
barrier has saturated to its maximum value /3 e g = 1/Tk- 
Instead fragile glasses have high excess complexity and 
still big variation of /3 e g along the supercooled line. The 
non-trivial temperature dependence of the activation bar- 
rier in (10) through the non-universal quantity (3 e s ex- 
plains why it is so difficult to find a unique empirical 
law that properly describes the viscosity anomaly of all 
glasses. The prediction that s* saturates to a finite value 
s* = Ap/Tx at Tk is not easy to check experimentally 
due to the difficulty to cover one order of magnitude 
in the effective barrier. Yet some experimental results 
suggest a crossover from fragile to strong behaviorEElEj. 
The present scenario could be also checked doing nu- 
merical aging experiments. According to ( |Io[ ) the t de- 
pendent effective temperature /3 e s(t) (measured through 
FDT-violations or the formula ||) should be given by 
&ff(t) -> (Tf/BA F )log(t/T ) if r a (t) ~ t and both r 
and BAp are nearly t and Tf independent. This rela- 
tion should hold for all Tf and t predicting a value for 
the activation barrier and the cooperative sizetia The 
experimental confirmation of a saturation of the hetero- 
geneity sizes for not too fragile glasses and the crossover 
from fragile to strong behavior is probably not out of 
reach and would be a check of the validity of the present 
theory. 
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